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PROBING THE EPOCH OF REIONIZATION WITH THE LYMAN ALPHA FOREST AT Z ~ 4 - 5 
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ABSTRACT 

The inhomogeneous cosmological reionization process leaves tangible imprints in the intergalactic 
medium down to z ~ 4 — 5. The Lyman-a forest flux power spectrum provides a potentially powerful 
probe of the epoch of reionization. With the existing SDSS I/II quasar sample we show that two 
cosmological reionization scenarios, one completing reionization at z — 6 and the other at z = 9, can 
be distinguished at ~ la level by utilizing Lya forest absorption spectra at z = 4.5 ± 0.5, in the 
absence of other physical processes that may also affect the Lya flux power spectrum. The redshift 
range z = 4 — 5 may provide the best window, because there is still enough transmitted flux and 
quasars to measure precise statistics of the flux fluctuations, and the IGM still retains a significant 
amount of memory of reionization. 

Subject headings: 



1. INTRODUCTION 

The history of cosmological reionization is presently 
primarily constrained by the cosmic microwave back- 
ground observation s of WMAP (W ilkinson Microwave 
Anisotropy Probe) (|Dunklev et al.ll2009D and the SDSS 
(Sloan Digital Sky Survey) quasar absorption spectra. 
The former gives an integral constraint, strongly sug- 
gesting that cosmological reionization may well be un- 
derway at z ~ 12, while the latter provides a solid anchor 
point at z ~ 6 when the universe beca me largely trans- 
parent to Lyman li mit photons (e.g. , IFan et al.l 120011: 
Becker et all 120011: ICen fc McDonald! 120021 : IFan eTatl 
2006bh . At z > 6.3 the lower bound on the neutral 
hydrogen fraction, x, of the IGM provided by SDSS 
observations is, however, fairly loose at x > 0.01. 
Thus, exactly when most of the neutral hydrogen be- 
came reionized is yet unknown and there are many 
possible scenarios that could meet the current observa- 
tional constraints (e.g. .iBarkana fc Loebl|2001l: ICeij|2003t 
Haiman fc Holder 2003: [Fan et al.ll2006at IWvithe fc Cenl 
20071 lBecke~e t al. 200J). 

The process of inhomogeneous cosmological reioniza- 
tion leaves quantifiable and si gnificant imprint s on the 
thermal evolution of the IGM (jTrac et al.ll2008f ). In this 
Letter, we show that the Lya forest flux spectrum at 
moderate redshift z — 4.5 ± 0.5 sensitively depends on 
and hence provides a very powerful probe of the epoch 
of reionization. 

2. REIONIZATION MODELS 

We use a hybrid code to accurately compute the reion- 
ization process, which consists of a high-resolution N- 
body code, a shock-capturing TVD hydro code and a ray- 
tracing radiative transfer (of Lyman l imit p hotons) code. 
The reader is referred to iTrac et all (|2008| ) for more de- 
tails. We use the best fit WMAP 5-year cosmological pa- 
rameters: il m = 0.28, fl A = 0.72, n b = 0.046, h = 0.70, 



o- 8 = 0.82, and n s = 0.96 (jKomatsu et al.ll2009l ). We use 
29 billion dark matter particles on an effective mesh with 
11, 520 3 cells in a comoving box of 100 h~ 1 Mpc, yielding 
a particle mass resolution of 2.68 x 10 6 h _1 MQ allow- 
ing us to resolve all atomic cooling dark matter halos. 
A total of N = 1536 3 gas cells of size 65kpc/h are used 
and we trace five frequency bins at > 13.6 eV with the 
ray-tracing code. The star formation rate is controlled 
by the halo formation history. We adjust the ionizing 
photon escape fraction to arrive at two models, where 
reionization is completed early {z ~ 9) and late [z ~ 6), 
respectively; note that the halo formation histories in the 
two models are identical. 

3. RESULTS 



Previous studi es (e.g., iFurlanetto et all 120041 : 
llliev et all l2006t iLee et alj |2008|) have shown that 
the reionization process proceeds in an inside-out 
fashion, where regions around high density peaks get 
reionized first. H II regions initially surround isolated 
galaxies that formed in high density peaks. With time 
these H II regions expand and lower density (void) 
regions are eventually engulfed by the expanding H 
II regions stemming from high density peaks. Conse- 
quently, the redshift of reionization of each individual 
spatial point, Zreion, is highly correlated with the 
underlying large-scale density field, with the positive 
correlation extending dow n to scales ~ 1 /i _1 Mpc, 
as we have shown earlier ([Trac et al.1 120081) . Once an 
expanding region is photo-ionized and photo-heated, it 
would cool subsequently due to adiabatic expansion and 
other cooling processes (primarily Compton cooling at 
high redshift), countered by photoheating of residual 
recombining hydrog en atoms (on the time scale of 
recom bination) (e.g.. iTheuns et aill2002l : iHui fc Haiman! 
2003!). As a result, the strong correlation between z re ion 
and the underlying large-scale density is manifested in a 
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Fig. 1. — Top panels show the log of the ratio of gas temperature from the simulation to that prescribed by a fixed EoS at z = 4, for 
the early (left) and late (right) reionization model, respectively. We use EoS formula T = 7b(p/po) ' 62 , where To is the temperature at 
mean density po m each model. The slice shown has a size (100 h -1 Mpc) 2 with a thickness equal to two hydro cells (130 /i _1 kpc). The 
distribution of flux transmission, F(early)= exp(— r(early)), for the late reionization model is shown in the bottom left. The flux difference 
between the two models: F (late)— F (early) = exp(— T(late)) — exp(— r(early)) is shown in the bottom right panel. 

strong anti-correlation between the temperature and the 
underlying large-scale density field. Specifically, different 
regions of the same low densities 8 < a few (without 
large-scale smoothing in this case) would display a 
large, long-range-correlated, dispersion in temperature, 
imme diately following the completion of reionization 
(e.g.. iTrac et al.l 120081 ). (Note that virialized regions 
are not affected and do not retain any information of 
reionization in this regard.) 

Both the anti-correlation between temperature and the 
underlying large-scale density and the consequent tem- 
perature dispersion at a fixed density weaken as time 
progresses and the temperature-density relation asymp- 
totically approaches a so-called equation-of-state (EoS), 
a one-to-one mapping from IGM density to temperature 



(|Hui fc Gnedinl[l99l . with T = T (p/p ) - 62 in the late- 
time limit. However, at the redshift range z = 4.5 ± 0.5, 
the IGM has not had enough time to have completely re- 
laxed to this state prescribed by the EoS such that quan- 
titatively significant deviations from a deterministic EoS 
exist, if the universe was reionized, say, at z r i ~ 6 — 8. 
The deviations from a simple temperature-density rela- 
tion are larger for smaller z r i at a given observed redshift. 

In Fig. [1] we show the log of the ratio of gas tem- 
perature from the simulation to that prescribed by the 
asymptotic EoS at z — 4 in a slice of size (100 ft, _1 Mpc) 2 
with a thickness equal to two hydro cells (130 ft. _1 kpc), 
for the early (top left panel) and late (top right panel) 
reionization model, respectively. The fields have been 
smoothed on cells of comoving length 130 /i -1 kpc. The 
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small reddish/yellowish regions seen in the top left panel 
correspond to virialized regions, for which the plotted 
ratio does not contain useful information. But these re- 
gions show clearly the location of ionizing sources. We 
see striking differences in temperature distributions be- 
tween the two reionization models with respect to their 
respective asymptotic EoS values. In the early reioniza- 
tion model (top left panel) most of the regions have blue 
color (i.e., the ratio equal to ~ 1) and appear to have 
mostly relaxed to the state predicted by the asymptotic 
EoS, while some low density regions in the voids still dis- 
play yellowish color with a temperature that is higher 
than that of the asymptotic EoS by 30 — 50%. On the 
other hand, in the late reionization simulation (top right 
panel), while regions just outside the shock-heated fila- 
ments and halos (bluish color) have largely relaxed to the 
asymptotic EoS, regions of comparable local densities in 
the voids are much hotter than that of the asymptotic 
EoS, by a factor of 1.5 - 2.5. 

Because the neutral hydrogen fraction in regions of 
moderate density is determined by the balance between 
photoionization rate and recombination rate, the latter 
of which is a function of temperature, the two different 
temperature distributions in the two reionization models 
result in different large-scale neutral hydrogen distribu- 
tion. In the bottom left panel of Fig [T]we show the ex- 
pected flux transmission, F(carly)= exp(— r(early)), for 
the early reionization model, where r(early) is the Lya 
optical depth computed based on the distribution of neu- 
tral hydrogen density, gas peculiar velocity and tempera- 
ture at z = 4 in the early reionization model. In comput- 
ing the neutral hydrogen fraction we have used a uniform 
background radiation field with its amplitude adjusted 
such that both models yield the same mean tr ansmitted 
flux of < F >= 0.43 at z = 4, as observed (|Fan et all 
2006b) . In the bottom right panel the flux difference be- 
tween the two models, F (late) — F (early), is shown, where 
it is clearly seen that the transmitted Lya flux is signif- 
icantly affected by the temperature difference at z = 4, 
resulting in fractional difference in the transmitted flux 
in the voids between the two models of ~ 15% (blue re- 
gions). Specifically, there is more transmitted flux in the 
void regions in the late reionization model, compensated 
by comparably reduced transmitted flux in high density 
regions. It is noted that, at z > 4, the majority of trans- 
mitted Lya flux comes from the lowest density regions 
of S < a few. 

Fig. [2] shows the ratio of flux power spectrum in the 
late reionization model to that in the early reionization 
model at z — 4 (black solid) and z = 5 (black dashed). 
It appears that the large-scale anti-correlation between 
density and deviations from a single EoS in the late reion- 
ization model leads to a significant amount of extra power 
in the flux spectrum (specifically, relatively high temper- 
atures in late-reionizing under-dense regions lead them to 
produce even less absorption than they otherwise would) . 
The difference between the flux power spectra of the two 
reionization models increases with scale, reaching 20% at 
k = 0.001(km/s) _1 at z = 4; the difference is still larger 
at z — 5 (~ 30%), as expected, due to still larger differ- 
ence in the temperature hence flux transmission between 
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Fig. 2. — Black solid and dashed curves are the ratio of flux 
power spectrum in the late reionization model to that in the early 
reionization model at z = 4 and z = 5, respectively. Also shown as 
the two green curves are the corresponding ratios produced by re- 
placing the real temperature in each simulation by that prescribed 
by the EoS given density (the same EoS in both simulations). The 
black error bars are the error one can expect from the full SDSS 
I/II sample plus existing high resolution data. The error bars will 
be approximately uncorrelated. A formal analysis of the 16 data 
points indicates that the two reionization models can be differen- 
tiated at 7<t level. 

the two reionization models. The black error bars in- 
dicate the statistical errors expected with the full SDSS 
I/II sample (completed, but not yet fully analyzed), With 
z ~ 4 SDSS I/II data plus existing high resolution data, 
one can distinguish formally between these two reion- 
ization models at 7a level. However, we note that the 
statistical differences between the two models are un- 
marginalized, i.e., not taking into account other physi- 
cal effects that a ffect the Lya flux pow er spectrum de- 
termination (e.g. iMcDonald et al.ll2005f h Therefore, the 
quoted statistical significance only serves as an indication 
of the potential power of this statistics. 

For comparison, the corresponding flux power spec- 
trum ratios at z = 4 and z = 5, if both models follow 
the same EoS (given density) , are shown (green curves) in 
Fig. [21 In this case, aside from the relatively small differ- 
ence on small scales due to cumulative dynamical affects 
on the gas density by the difference in the gas pressure 
histories, the two models have identical flux power spec- 
trum on large scales. This clearly demonstrates that the 
large difference in the flux power spectra between the two 
reionization models (black curves in Fig. [5J is a result 
of large differences in the contemporaneous temperature 
distributions. 

4. DISCUSSION 

This effect of inhomogeneous reioniz ation on the flux 
power spectrum was explored earlier bv lLai et alj (|2006f ) 
at z = 3, based on a semi-analytic model. Their fo- 
cus is on z = 3 and found that, on large scales, k ~ 
0.001(km/s) , temperature fluctuations lead to an in- 
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crease in the z ~ 3 flux power spectrum by at most 10%. 
Our focus here is at higher redshifts z = 4 — 5 and the 
effects, not surprisingly, are larger and potentially more 
discriminating. 

A fluctuating radiation background, produced largely 
by radiation from sparsely distributed quasars but 
also by galaxies, can a ff ect the flux power spectrum 
(iMeiksin fe White! l200i ICroftl [200l IMcDonald et~aTl 
2005). Larger fluctuations in the radiation background 
give rise to larger a mplitudes of the flux pow er spectrum 
at large scales (e.g. IMcDonald et al.l 120051 Figures 6,7 
therein). This enhancement of the flux power spectrum 
on large scales due to a fluctuating radiation background 
will be in addition to what is caused by the gas tem- 
perature fluctuations shown here, if QSOs were dom- 
inant. The radiation contribution from stars may be 
more dominant at the redshi ft range of concern here (e.g., 
iFaucher-Giguere et al.l l2009) . Star formation is known to 
be biased and hence higher density regions, on average, 
tend to have higher radiation field than lower density re- 
gions. Thus, the two effects due to a fluctuating radiation 
background and an inhomogcncous reionization process 
may be partially degenerate or have a tendency to cancel 
each other's contribution, although there is a possibility 
that the radiation fluctuations may b e relatively mod- 
est (e.g., Mesi nger fc Furlanettol [2009( 1 . A more careful 
modeling of the contribution from quasars as well as radi- 
ation sinks (such as Lyman limit systems) is required in 
a comprehensive modeling. The purpose of this Letter is 
to demonstrate that, if the effects on the Lya flux power 
spectrum determination due to the epoch of reionization 
were the only relevant ones, then a precise measure of 
the flux power spectrum with the full SDSS 1/11 data 
will be able to place a very tight constraint on the epoch 
of reionization. 

However, a detailed comparison between models and 



SDSS I /ll observations requires a full analysis of all 
astrophysical/cosmological processes that may affect 
the determination of the flux power spectrum and 
some of them may be degenerate to varying degrees 
([McDonald et al.l l2005L including fluctuating radiation 
field, damped Lya systems, galaxy formation feedback, 
initial photoheating temperature (i.e., related to IMF of 
high redshift galaxies), X-ray heating, He II reionization, 
among others, before its statistical potential can be pre- 
cisely marginalized and quantified. We will perform such 
an analysis in a future study. 



5. CONCLUSIONS 

Utilizing state-of-the-art radiative transfer hydrody- 
namic simulations of cosmological reionization, we put 
forth the point that the inhomogeneous reionization pro- 
cess imprints important and quantitatively significant 
signatures in the intergalactic medium at z = 4.5 ± 0.5 
that can be probed by the Lya forest in the quasar ab- 
sorption spectra. We illustrate that with Lya forest data 
at z = 4 - 5 to be provided by the SDSS I/II full data 
sample, one may be able to distinguish between two cos- 
mological epochs of reionization, one at z = 6 and the 
other z = 9 at 7cr level, if they were the only effects on 
the determination of the Lya flux power spectrum. 
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